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Abstract. We obtain LP eigenfunction bounds for the harmonic oscillator 
H = — A+x 2 in IR ra and for other related operators, improving earlier results 
of Thangavclu and Karadzhov. We also construct suitable counterexamples 
which show that our estimates are sharp. 



1. Introduction 

The question of obtaining LP eigenfunction bounds for elliptic operators on 
compact manifolds has been considered in Sogge's work, for which we refer the 
reader to his book [TT]. The LP eigenfunction bounds in are sharp, and turn 
out to be related to a variable coefficient version of the restriction theorem, 
and further to a phase curvature condition for Fourier integral operators. In 
this analysis a special role is played by the Laplace-Beltrami operator for 
the sphere, which is the worst case because it has many highly concentrated 
eigenfunctions. This is connected to the fact that it has a periodic Hamilton 
flow. 

In this article we consider the problem of obtaining LP eigenfunction bounds 
for the Hermite operator LL = —A + x 2 in M. n and also for a larger class of 
related operators of the form Hy = —A + V. Within this class the Hermite 
operator plays a role similar to that of the spherical Laplacian, in that it has 
a periodic Hamilton flow and many highly concentrated eigenfunctions. 

This question has received considerable interest in the context of Riesz 
summability for the harmonic oscillator in the work of Thangavelu ^2], |14j . 
[T3] and Karadzhov jjj. 

Our interest in it has a different source, namely the strong unique contin- 
uation problem for parabolic equations. In this context the work of Thangavelu 
and Karadzhov has already found applications in Escauriaza j2j and Escauriaza- 
Vega Further applications are contained in a forthcoming paper of the 
authors. We note in passing that the related strong unique continuation prob- 
lem for second order elliptic operators is related to the eigenfunction bounds 
for the spherical harmonics. This was first observed in work of Jerison [0] ; see 
also the authors paper [TO] and further references therein. 
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In the next section we begin by considering the problem of obtaining dis- 
persive and Strichartz estimates for the corresponding Schrodinger equation. 
This leads to an alternative proof of the eigenfunction bounds of Karadzhov j^j 
and Thangavelu ^Hj- Our approach has the advantage that is robust enough 
so that it allows us to obtain the same results with x 2 replaced by potentials 
in a very large class. 

Then we direct our attention to the eigenfunctions (— A+x 2 )4> = A 2 0. These 
are concentrated inside the ball {|x| < A}, and have an exponential Airy type 
decay beyond this threshold. The behavior of the eigenfunctions inside the 
ball is not very different from (a rescaling of) what happens in a bounded 
domain. However, considerable care is required near the boundary of the ball, 
where the concentration scales are different. Consequently, it is more natural 
to obtain weighted LP estimates with weights which are essentially powers of 
A — \x\. The results we obtain strengthen those of Karadzhov and Thangavelu 
and complete the picture. As before, our methods are robust and apply equally 
to any potential which behaves roughly like x 2 . 

In the last section of the paper we construct appropriate examples which 
illustrate the possible concentration scales for eigenvalues of the Hermite op- 
erator and show that our LP bounds are sharp. 

To conclude the introduction we provide the reader with a special case of 
our main result. Denote the spectral projection to the eigenvalue k = A 2 by 
Pk- Then our LP bounds for eigenfunctions of the Hermite operator H have 
the form 
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2. Dispersive estimates for the Schrodinger equation 
Consider the Schrodinger operator associated to the Hermite operator in 

id t -H 

which generates a group of isometries t — > e ttH . Furthermore, this group is 
periodic with period n. Here we investigate the question of obtaining pointwise 
bounds for the kernel of e ltH . We also consider the same question for more 
general operators 

H v = -A + V 
Theorem 1. a) The operators e ltH satisfy 

||e itH ||L^£« < |sint|-5 

b) Let V be a potential which satisfies 

(2) \d a V\ < c a , \a\ > 2 
Then for \t\ <C 1 the operators e ltH satisfy 

\\ itH v \\ < 

Proof of TheoremU\ Prove the bound at time to- 111 ca se (a) the estimate 
follows from well known formulas for the Schrodinger kernel of the Hermite 
operator. We prefer to give a more flexible argument. By periodicity we can 
assume that |to| < tt/2 and replace | sin 1 1 by |£|. We rescale to time 1 by 
setting 

t x 

t ^ ^ X ^ ~. — 

^0 v^o 

Then we need to prove an uniform bound for the kernel of e lHt ° , respectively 
e tH v where 

H t0 = tlx 2 - A, Hp = t y(iA) - A 
Since to is i n a bounded set, it follows that the symbols 

h to (x, o = t 2 x 2 + e, h#(x, o = t v(xvr ) + e 

satisfy the bounds 

(3) \d^dlh{x,i)\<c a ^ \a\ + \(3\>2 

uniformly in |io| % 1- This implies that we can use directly Proposition 4.3 in 
[H] to obtain a phase space representation of the fundamental solution K(t, y, y) 
for idt — Hy. For the reader's convenience we restate the result here: 
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Lemma 2.1 ([H])- Let h be a symbol which satisfies TTien for \t\ < 1 i/ie 
fundamental solution K(t, y, y) for id t — h w (x, D) can be represented as 

(4) K(t,y,y)= [ e-^'^e^-^e^^e^^G^x^y)^^ 



where (x, £) to t/ie Hamilton flow for h, the function G satisfies 

(5) | (a* - yy%dl%G(t 7 x,t,y) I < 

and i/ie rea/ phase function if) is determined by 

= -Kx^e) +eh(x t ,e), = o 

A feature of the construction in jH] is that the integrand solves the evolution 
equation for each (x,£). In addition, it is concentrated in the phase space on 
the unit scale along the bicharacteristic t —>■ (x t ,^, t ). Such highly localized 
solutions are called wave packets. With this terminology, one can view the 
above lemma as a way of representing solutions for id t — h w (x,D) as almost 
orthogonal superpositions of wave packets. 

In our case we need a pointwise bound for K(l, y, y). Neglecting all oscilla- 
tions in (j3J) we write 

\K(l,y,y)\< [ e-^-^l + ly-x 1 !)-^^ 



Due to (JHJ) the Hamilton flow for Hy is Lipschitz. Hence the integration in x 
is trivial, and we obtain 



\K(l,y,y)\< [ + 



In order to obtain an uniform bound it suffices to prove that the Lipschitz map 
£ — ► y 1 has a Lipschitz inverse. In the case of the Hermite operator this map 
is linear, and the desired conclusion is obtained by direct computation. For a 
more general potential V one computes the linearization of the Hamilton flow, 
which shows that 

^L = 2I n + O(t ) 

For small to this shows that the map £ — >• y 1 is a global diffeomorphism of 
R n . □ 

As a consequence of the dispersive estimates for the Schrodinger equation 
we also obtain Strichartz estimates for the Schrodinger equation: 

Theorem 2. Let V be a potential which satisfies 
(6) \d a V\ < 1, \a\ = 2 

Then the solution u to 

{id t ~ H v )u = f, u(0) = u 
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satisfies 

\\u\\ L pi(o,i;m) < \\u \\ L 2 + ll/ll^^.^) 
whenever the pairs (pi,qi) and (^2,^2) o^e subject to 

2 TL — 1 TL — 1 

-H = — =— , 2<p<oo, (n,p,g) ^ (2,2,oo) 

p q 2 

Proof, a) If V satisfies the stronger bound (J2J) then this follows from Theorem^ 
by standard arguments as in Ginibre-Velo [S] and the references therein. For 
the endpoint one can use the results in Keel-Tao jH]. 

b) For potentials V which only satisfy (jUJ) we use a frequency decomposition 
of V. Given a smooth compactly supported function x which equals 1 in a 
neighborhood of the origin we set 

V = V + V u V Q = X (D)V, Vx = (1 - x(D)V) 

The low frequency part V satisfies J2J), therefore we use part (a) of the 
proof. On the other hand the high frequency part V\ is bounded, so we can 
add it in using the fact that the result in the theorem is stable with respect 
to L 2 bounded perturbations of Hy. □ 

Using the Strichartz estimates one can easily obtain eigenfunction bounds. 
We begin with the Hermite operator H . Its spectrum is n + 2N. A basis of 
orthogonal eigenfunction is given by the Hermite functions, 

K = (-l)\<*\ e x2 l 2 d a e- x2 

with 

Hh a — (n + 2|a|)/z Q 

For an integer k = X 2 in the spectrum we denote by Pk the corresponding 
projector. For G L 2 the function 

u = e uk P k <J) 

solves (idt — H)u = 0. Hence we can apply Theorem |2 to obtain: 

Corollary 2.2. Let P^ be a spectral projector for H . Then 

2n 

WWlp < U\\l*, 2<p< ——, (n,p) ^ (2, 00) 



In particular this gives L p bounds for each eigenfunction. The estimate still 
holds in the case (n,p) = (2, 00), but to prove this one needs the stronger 
results in the next section. 

In the case of the operator Hy the spectrum need not be discrete. Even if it 
is, there is no guarantee that there is an 0(1) spectral gap. Hence it is more 
useful to consider spectral projectors P[k,k+i] associated to unit size intervals. 
For <fi G L 2 the function 

u = e itk P [Kk+1] ct) 
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satisfies 

\\(id t - H)u\\ L <*>& < ||0|| L 2 
Hence we can apply Theorem El to obtain: 

Corollary 2.3. Let V be a potential which satisfies JBJ). If P\k,k+i] is a spectral 
projector for Hy then 

2n 

\\P[k,k+i]<l>\\L P < W\\v, 2<p< -—^ (n,p) ^ (2, oo) 

Here and in the sequel we formulate results for n > 2, which are valid for 
n — 1 with modifications which are either obvious or discussed. 

3. Weighted W eigenfunction bounds 

Let (0, A 2 ) be an eigenfunction, respectively an eigenvalue for the Hermite 
operator. The function is essentially concentrated in the ball {\x\ < A} 
(modulo an exponentially decaying tail). We split the interior of this region 
into overlapping dyadic parts with respect to the distance to the boundary, 

Df = {\x\e [A(l - 2- 2 ^ 1 )), A(l - 2~^' +1 >)]}, 1 < 2 j < A§ 

By D bd we denote a narrow strip near the boundary of the ball, 

D bd = {\\x\-X\ < A~s} 

Finally, D ext is the exterior of the ball, 

D ext = {\x\ > A + ^A-5} 

The symbol of x 2 — A — A 2 equals 

x 2 + e 2 -A 2 

If x G DJ 1 * then this symbol can be zero only in the region 

ICI « A2"' 

2 1 

Hence the case 2 ] « As corresponds to an annulus of thickness A - 3 and to 
|£| rti As which is exactly on the scale of the uncertainty principle. This 
explains the limitation for the range of 2 3 . It also gives the scale on which (f> 
decays away from the ball. 

Given an eigenvalue A we define the spaces l q x L p of functions in R n with 
norm 

l<i,2^<A 2 /3 3 

with the usual modification when q = oo. The subscript A is used because the 
sets DJ 1 *, D bd and D ext depend on A. 
For i 6 K" we denote 

y = A~§(A 2 -x 2 ), (y)_=l + y_, (y)+ = l + y+. 
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Then 

Theorem 3. a) Let 2 < p < Then 

(7) ||A^t^-)( y )7 + 4 ^~'%^<U\\v' + \\(H-X^\\v> 
b) Let < P < oo. T/ien 

-, 11 i H f i _ 2_ 'l 

(8) IIA'-^-^r* (2/>-^0||^ < 

In the interesting regions Dj the estimates can be reformulated as follows. 
If 2 < p < then 

— r — n—l 

(9) xt-nfr-^G-'HttlwDF) < M\\* + W - \*)<p\\ L > 

and if 2 (" +1 ) < p < oo then 

n— 1 — * — 

(10) X 1 ~ n ^~p ) 2 j{1 ~ n ^~p^\\P X 2(f)\\ LP{D ^t) < ||0|| L2 

The estimates in Theorem El are sharp. We give here a heuristic motivation, 
and complement this with a more precise (yet not as general) construction in 
the last section of this paper. 

The symbol of H is x 2 + £ 2 . The A 2 eigenfunctions of H are concentrated 
in the phase space within a neighbourhood of size 1 of the characteristic set 
{x 2 + £ 2 = A 2 }. For a more precise description, one can consider the Hamil- 
ton flow for H restricted to the same sphere. It is periodic, and one can can 
construct eigenfunctions which are concentrated in the phase space within a 
neighbourhood of size 1 of each such bicharacteristic. To obtain eigenfunctions 
which are pointwise larger one may consider superpositions of such eigenfunc- 
tions corresponding to neighbouring bicharacteristics. However, such a con- 
centration can only occur on a smaller set in the physical space; the optimal 
balance between the amplitude and the localization region is dictated by the 
uncertainty principle. 

Consider first the estimate J7J). Within D™ 1 the worst concentration occurs 
for eigenvalues concentrated both in space and in frequency in a 1 neighbor- 
hood of a single bicharacteristic. For an arbitrary j one needs to consider 
modes which within DJ 1 * are spatially concentrated within 2~2 of a bichar- 
acteristic, with a dual concentration in frequency. In the extreme case, for 

2 1 

2 J = As, the concentration occurs as a bump in a A~3 ball near the circle 
\x\ = A. Except in dimension 1, this corresponds to a solution which is uni- 
formly spread in an ellipsoid of size A x As. 

A different behavior is responsible for estimate (jSJ). Here the concentra- 
tion corresponds to an equidistribution of energy between all bicharacteristics 
through a point P. For P G DJ 11 this leads to concentration in a 2 J A~ 1 ball 
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around the point. If for instance P is the origin then this could be a spherically 
symmetric mode. 



Remark 3.1. We expect that the estimate holds with l^L p replaced by LP for 
all p > 2. This is because the concentration cannot occur simultaneously in 
many regions DJ 11 . From Corollaru \2.^\ we already know that the improvement 
holds if ^2 — P — 00 ■ On the other hand, such an improvement cannot be 
true for p = 2. 

If we want to write an estimate without the weights then the sign of the 
power of (y) + becomes important: 

Corollary 3.2. Let (<f), A 2 ) be an eigenf unction, respectively an eigenvalue for 
the Hermite operator. Then 

where for n > 2 we have (see Figure^) 
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In the first and the third case the worst bound is in D™*, while in the second 
case the worst bound is within D bd . 

Endpoints discussion: The arguments below yield the p = oo endpoint even 
in dimension n = 1,2. For n > 3 the case p 



-^r comes from the Strichartz 

n—2 



estimates in the previous section. The p 



2(n+3) 
n+1 



endpoint is false for n 



where the eigenvalues are simple and have an Airy type behavior at the ends. 
However, it is likely be true for n > 2, because it is not possible to have 
concentration at all scales. 

Proof of Theorem^- The L? bound. This is the key step in the proof of The- 
orem El because the L? bound is strong enough to provide the localization 
which is needed for the rest of the arguments. We assume that 



l 2 + ||(#-A< 
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and we shall prove the following bounds: 

U\\l\d^) < 2"5, ||V0|| L2(jD « t) < A2"^ 

(11) U\\lhd^)<^K ||V0|| lW <1 

\\{x 2 - A 2 )0|| L2( ^) < A*, ||(x 2 - A 2 )^V0|| L2(D -) < A^ 

We first offer some intuitive justification for these bounds. The one in D ext 
is an elliptic estimate. For the rest, it is convenient to think of an almost 
orthogonal basis for the eigenspaces of H which consists of eigenfunctions 
which are localized in the phase space within a neighbourhood of size 1 of null 
bicharacteristics of H — A 2 . The energy of such eigenvalues is equidistributed 
along the corresponding bicharacteristics, which have length A. Hence in order 
to measure what portion of the energy is contained in DJ 1 * it suffices to measure 
the length of the bicharacteristic within DJ 1 *, which is at most 2~ J A. This 
justifies the L 2 bounds. For the gradient bounds we simply note that on the 
characteristic set of H — A 2 within Dj we have |£| ~ 2~ J A. 

To prove we begin with some simpler estimates. First note that 

(0, (i/-A 2 )0) = ||V0||i 2 + ((x 2 -A 2 )0,0) 

Since x 2 — A 2 > —A 2 , this easily leads to 

(12) l|V0|U 2 <A 
This argument can be easily improved. Choose 

W = ((\ 2 - x 2 ) 2 + Xl)- 1 * 
Then an integration by parts yields 

{W<j>, (H - A 2 )0) = \\W^m 2 + 0(U\\ L2 ) 
which produces another nonsharp estimate, 

(13) ll((A 2 -x 2 ) 2 + At)-3V0|| L2 <l 

which is slightly better than |T2|) near the sphere \x\ = A. We use this weaker 
bound to eliminate some error terms in the estimates which follow. 

We will argue in a manner which is similar to the Carleman estimates. 
First we introduce a bounded exponential weight which does not change the 
estimates, but allows us to replace the operator H — A 2 with its conjugate with 
respect to the weight, namely the operator H a ^\ below. The advantage is that 
the conjugate operator is no longer selfadjoint; precisely, the gain in the L 2 
estimates comes from the positivity of the commutator between its self-adjoint 
and its skew-adjoint part. In order to guarantee this the weight needs to be 
chosen roughly so that it is convex along the null bicharacteristics for H — A 2 
near the ball {|x| = A}. 

An alternate approach would be to obtain a Morawetz type estimate using 
a suitably chosen multiplier. We do not pursue this as it seems slightly less 
precise, and requires more care in the error estimates. 
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Consider a weight function a : R — > K. with the following properties: 

(i) a is nondecreasing and equals in (—00, —2]. 

(ii) a and its derivatives satisfy the bounds 

(14) \d k a(y)\ <(l + | y |)^- fc , k>0 

Let c be a small positive constant which will be chosen later. All implicit 
constants in the estimates that follow will be independent of c. Denote 

j,(x) = e- x ~ ia ( y U(x), V = c\-l{\ 2 - x 2 ) 
Introduce also the conjugated operator 

H a \ = e- x ' ia{y \H - A 2 )e A ~ §a(y) 

Then 

H a ^ = e- x ^ a ^(H-\ 2 ) ( p 

The exponential weight is bounded since x > implies y < A3. Hence 

U\\ L 2 + \\(H- A 2 )0|| L2 « ||^|| L2 + \\H a M» 

On the other hand, given the above properties of a it is easy to verify that ip 
satisfies fill)) if does. Hence we have replaced <fi and H — X 2 by ip, respectively 
H a ,\- 

The conjugated operator H a ^\ is decomposed into a selfadjoint and a skew- 
adjoint part, 

tt z_rre , Trim 

"a, A - -n a ,A + ^a,A 

iZ™ = -A + x 2 - A 2 - 4c 2 A-ix 2 |a'(y)| 2 

flS = -2cA-f(9a;a , (j/) + a , (y)x9) 

where we use a short notation for operators dxa'<j) = Y^i=i 9 Xi {xiCi'4>)- 

Ideally we would like the commutator {H™ X ,H™} to be positive definite. 
However, this is too much to hope for. Instead, it turns out that we can 
choose the exponential weight so that the commutator has a positive symbol 
only near the characteristic set of H a) \. To compensate for that we introduce 
a real correction term W\ of the form 

W x = 4c 2 b(c\-I(\ 2 -x 2 )) 

where the positive function b is chosen to satisfy the bounds 

(15) \d k b{y)\<(y)-(y)7(y)-\ k > 

Then we compute 

\\H a ^\\ 2 L2 = ||^||i 2 + ||^>||i 2 + (K, e A ,^]^^) 
= \\(h% - WxMh + WK^Wh + ^ 

where 

c = [h:% h™\ + h:%w x + w x h:% - w 2 
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Therefore we obtain 



(16) (Cij,i')<\\HaMh<l 

Our goal is to choose the weights a and b so that C is a positive operator which 
controls the norms in for ifj. Hence we need to compute C. This can be 
rather tedious, and in order to simplify the analysis we make two conventions: 

(a) We discard all derivative terms which can be controlled by (J13|) . 

(b) All the scalar terms which are negligible are incorporated into a generic 
term called "error", which satisfies 

|error| < 1 + (?\%{y)+ l (y} 2 _ 

As we shall see later on, the error term in (b) is easily controlled by the 
main term in (|17|) below. For instance the scalar W 2 is an error term. We 
consider the two remaining terms in C. Using the bounds on the derivatives 
of a in (ii) above yields 

[H r a % H l a m x ] = 8c{\-lx 2 a'(y) - 2c\~ 2 dxa" \y)xd + \-*da'{y)d) + error 
On the other hand 

W x H r a % + H r a %W x = -8c 2 (db(y)d + (x 2 - A 2 )%)) + error 
Summing up the two sets of estimates and using Cauchy-Schwartz we obtain 

1 > j ($c 2 b{y) - 16c 2 A- 2 xV (y) - 8c\-*a'(yj) |V0| 2 

+ (8c\~^x 2 a'(y) - 8c 2 b(y)(x 2 - A 2 ) + error j \(p\ 2 dx 

where a"_ is the negative part of a", a"_ = (\a"\ — a")/ 2. 

Using the estimates (|12j) and (|13j) we can control the third derivative term, 
and also replace the x 2 in the second derivative term by A 2 . The x 2 in the 
first scalar term can be replaced by A 2 modulo a bounded remainder. Thus 
we obtain 

(17) 1> j (c 2 (b(y)-2a'Uy)))\V<j ) \ 2 +(c\l(a'(y)-yb(y))+ error) \<p\ 2 dx 

To conclude the proof we need to choose the functions a and b so that the 
two coefficients above are nonnegative and sufficiently large. The key step is 
summarized in the next Lemma. 

Lemma 3.3. Let {Ek)k>i be a slowly varying positive sequence with the prop- 
erty that Yli £ h — 1- Then there are functions a, b as above so that in addition 
we have: 

{ (i + M) y<i 

b(y) - 2a"_{y) > 

( e k \y\~l ye [2\2 k+l ] 
ii 



a'(y) - yb{y) ^ { 

( e fc |y|-a y e [2\2 k+l ] 

We first show how to conclude the proof of the L 2 estimates (jllj) . First 
we need to control the error term, and this is where we use the freedom to 
choose c sufficiently small. This is because the positive scalar term in (I17J) has 
a factor of c, while the unbounded part of the error has a factor of c 2 . Hence 
in order to control the error it suffices to insure that 

- k 

e k > 2~ 

which guarantees that 

a'(y)-yb(y) > 

The bounds in the lemma suffice to obtain in the regions D ext and D bd . 
In order to prove (fTTj) in the region DJ 1 * we need to also choose the e^s so 
that 

£fe ~ 1, when 2 k+2j « As 

Proof of Lemma \3.tft (i) The range y < 0. Observe first that we can add any 
positive function to b for y < and improve both inequalities. Hence as long 
as a'(0) > and a"(0) > we simply need to choose b sufficiently large and 
growing like \y\ at — oo. 

(ii) The range y > 0. Here it suffices to choose a so that a' > 0, a"{y) < 
for large y, a"(0) > and 

a\y)-2ya"{y)>e k \y\^ ye[2\2 k+1 ] 

Then b can be chosen b(y) = a"(y) + (2y)~ 1 a'(y) for large y and arbitrary 
between 2a"_{y) and y a'(y) for small y. 

It remains to describe the choice of a'. We begin with a simple choice, 
namely 

a' (y) = (1 + y 2 )"i 
which satisfies all conditions except for a weaker bound 

a' (y) - 2ya>>(y) > (I + y 2 )^ 

Then we choose a nonincreasing function <i with | < d < | so that cf(0) = 
and 

<f(y) ~ -efci/" 1 y G [2 fc , 2 fc+1 ] 

Finally, we set 

a'(y) = d{y)a' Q {y) 

Then in [2 k , 2 k+1 ] we have 

a'{y) - 2ya"(y) = d(y)(a Q (y) - 2ya'o(y)) - 2ya' (y)d'(y) > e k \y\~^ 

□ 
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Proof of Theorem 0' the L bound. We split -D* n * into balls £)* n ' ,fc Q f 

radius Tj = and consider a corresponding partition of unity 

i = x + x + Xj ( x ) 

with 

\^X%x)\<rj\ \V\%x)\<rf 
The localized pieces of (f) are 

$ = X -<P, <P hd = x b % <P ext = x ext <P 

We consider three cases. 

a) The interior estimate. We claim that the functions (f>^ satisfy 

(is) E^n^ii^+^ii^-^iii' 2 ^ 1 

k 

The first part follows directly from the lemma, for the second we commute 

(H - A 2 )0j = [-A, x*]<f> = V^V0 + V 2 X ?0 

and use the lemma and the bounds on Xj- 

To obtain LP bounds for each of these pieces we use a Strichartz type esti- 
mate which is a special case of Theorem 2, jHj. For convenience we state it in 
the following 

Lemma 3.4. Let W be a real potential in the unit ball which satisfies W ~ 1 
and 

\d a W\<l, |a| = l,2 

and let a* J be elliptic coefficients of class C 2 . Then for all \i > 1 and u 
supported in the unit ball we have 

\\u\\ 2(n+i) < n~^{^\\u\\ L 2 + fi'^lKa^didj + fi 2 W)u\\ L 2) 

L 



Discussion. Strictly speaking, in order to apply Theorem 2, |U] we need the 
additional bounds 

\d a W\<ii 2 ¥ L , |«|>2 
But such bounds can be easily gained by truncating the potential W in fre- 
quency, 

W = W + W 1 W = x{^D)W, W x = (1 - x{^D))W 

Then the low frequency part Wq satisfies the stronger bounds. The high fre- 
quency part on the other hand satisfies a pointwise bound \W\\ < /i" 1 there- 
fore it does not affect the size of the right hand side. Another observation is 

2(n+l) 

that Theorem 2, [2] only gives the L bound for u in the frequency re- 
gion |£| < fi. However, outside this region the symbol — £ 2 + [i 2 W is elliptic, 
therefore even stronger bounds are easy to obtain. □ 
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r, _____ 2(n + l) 

Given the L bound (jI8|) . the L estimate for in DJ 1 would follow 
from an estimate for (p 1 *, namely 

U-W ^±11 <(2^ 1 )^(2^||^|| L2 +2^||(-A + x 2 -A 2 )^|| L2 ) 

£ 71 — 1 J J 

This follows from Lemma 13.41 after rescaling to the unit spatial scale. Note 
that within £)*. nt ' fc the symbol x 2 + £ 2 — A 2 is elliptic at frequencies |£| ^> 2 -J A, 
so the interesting region in phase space has size (2 _2j A) n x (2~^X) n . After 
rescaling, the frequency becomes \i = 2~^\ 2 > 1. 

b) The boundary estimate. We can use the bounds in to obtain 

Ai||0 M || i2 + ||V0 M |U- <1 

2(n + l) 

Then the L bound for is straightforward by Sobolev embeddings. 

c) The exterior estimate. The bound in (fTTj) implies that 

|| ( x 2 _ x 2 )<j> ext \\ L , + || (x 2 - A 2 )^V<r1 L2 < 
Then by (weighted) Sobolev embeddings we obtain 



||(x 2 -A 2 )w>V<ri < A3 

□ 

Proof of Theorem^- the L°° bound. We consider separately two cases. 



(i) In D mt and D bd we only use the size of the potential. Each region D v - 

has size rj = 2~ 2 i\ and corresponds to frequencies of size jij = 2~i\. In D % - 
we have an elliptic estimate, 

n{~--) 2p 

" P llxj^llw ^ ^7 2 |l A 0ll L p(£.f*' fe ) + UW^Df^y 1 < P < g < 



int,k 



int,k 



n — 2p 

which follows by rescaling from the case r = 1, \i > 1. 

We replace A0 with (A — x 2 + A 2 )0 + (x 2 — A 2 )0, begin with p = and 
apply the above estimate iteratively until we arrive at q = oo. The same idea 
works also in D bd . 

(ii) In D ext we also use the favorable sign of the potential. We first improve 
the L? bound on 0, namely the last part of We rewrite (jllj) as 

ll(2/)+0lU 2 (D-') < A~^ 

and inductively show that 

a i 
(19) ||(y)|0|| i2 pe, t) <A-3 

Suppose that (JTHJ) holds. Then we compute 

x ext (y)T<f> (y)J(H - \ 2 )<l>dx = J x ext {y) N + {{x 2 - a 2 )</> 2 + |V0| 2 ) 

+ x ex "(i/)^- 1 A-^xV0 + w x ext (y)^<t>dx 
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In the region y ^> 1 the last two right hand side terms are controlled by the 
first. In the region y ~ 1 we use (jlip . Together with (J19j) this yields 

J X e ^(^((x 2 -A 2 )0 2 + |V0| 2 )rfx< 1 

which implies (II 9|) with N replaced by A^ + 1. 

At this point we can conclude the proof as in case (i). Precisely, to fix the 
size of the potential we also consider a covering with balls 

jjext ^- jjext,k 
j>l,k 

and a corresponding partition of unit 

1 = X) in D ext 

so that in D^ xt,k we have 

\x\ 2 - A 2 ~ 2 2 ^A§ 

In each of these balls we use the same elliptic argument as in (i). □ 

4. Extensions 

In the previous section we used the potential x 2 . However, its precise form 
does not play a fundamental role in the estimates. Here we consider instead 
the operator Hy with a positive potential V which satisfies the following con- 
ditions: 

(20) V~\x\ 2 , [W| ~ |ar[, \dlV\ < 1 

Given an eigenvalue A 2 of V we introduce as before the dyadic regions 

D int = | A 2 _ y e p-20-l) A 2 ; 2 -20'+ 1 )A 2 ] 2 < 2 j < A§} 

D bd = {\\ 2 -V\ < A§} D ext = {V > A 2 + A§} 

Also we set 

y = \-i( V -\ 2 ) 

Then 

Theorem 4. a) Let 2 < p < Then 

(21) ||A^t( 5 -) (y) 7 + * (5 "' J (y> + " 7t3 "'Vlli^^ IHb + ||(^-A 2 )0|| L2 
frj Lei 2 ^ t J f 1 1 - > < p < oo. T/ien 

, , , 1 n I 1 _ 1_\ 

(22) llA^a^-^^r^^^^^^+H^II^ £ IWI^ 
The result in Corollary 13.21 also applies to Hy. 
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Proof. (Sketch) We only need the counterpart of (JTTJ), the rest of the proof is 
identical. 

Without any restriction in generality we can replace V with a mollified 
potential x(D)V, since the difference is bounded. This allows us to use the 
additional assumptions (J2J). 

We use the same functions a and b as in the proof of (fTTj) but with the above 
definition of y. Now the selfadjoint and skew-adjoint parts of the conjugated 
operator are 

H r a % = -A + V - A 2 - c 2 A-!|W|V(y)| 2 

H l a m x = -c\-HdVVa\y) + a'{y)VVd) 
The correction term is chosen to be 

W x = 2c 2 A" 2 |W| 2 %) 
Modulo error terms the operator C has the form 

C = 2c(\-l\VV\ 2 a'(y) - A~ 2 %)| W| 2 (V - A 2 )) + 2c 2 d(\- 2 \VV\ 2 b(y) 
+ \- 2 VVa"{y)VV)d 
Then the argument continues as in the proof of (fTTj) . 

□ 



5. Remarks on optimality 

We begin by recalling bounds and expansions of Hermite functions. The 
Hermite functions are the eigenfunctions of the one dimensional Hermite op- 
erator, and solve 



(23) -h'l + x 2 h k = {2k + l)h 

for nonnegative integers k. They are given by 



dx K 



and are even functions for even k and odd functions for odd k. Here the 
meaning of k differs slightly from the previous sections. In dimension n a 
complete set of eigenfunctions is given by 

n 

h a (x) = Y\_ h a t (Xi) 

i=l 

where the corresponding eigenvalue is n + 2\a\. 

To construct highly localized eigenfunctions we need a better understanding 
of the behavior of the Hermite functions. This is well understood by now, and 
we describe it next. 

The ODE has a turning point at x = \/2k + 1. We set A = \/2k + 1. 
Then the functions hf. have an oscillatory behavior for small x, an Airy type 
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behavior for \x\ close to A and Gaussian decay for large x. More precisely, 
define 

(24) s~(x) 



Jo 

J* Vl* 2 -A 2 |dt, 



(25) s + (x 
Then 

Lemma 5.1. The normalized eigenfunctions 

hk = hk\\hk\\ L 2 

satisfy 

°2fc(A 2 — x 2 )~i (cos s~(x) + error) \x\ < A — A - ^ 

/^2fc = <( 0(As) A- \~\ < x < A + A _ 3 

a^ fc e~ s+ ^(A 2 — x 2 )^3(l + error) \x\ > A + A~s 

a 2fc+i(A 2 — x 2 )~3(sin s _ (x) + error) \x\ < A — A~s 



^2fc+i = "J 0(A«) A - A"3 < x < A + A~^ 

l 2fc+l 



3 S+ ( :C )(A 2 — x 2 ) 4(1 + error) x>\ + \ 3 



where 



,±1 



1, error = 0((|x 2 — A 2 | 2 ||x| — A| x ) 

Note that the error term is 0(1) if — A| ~ A~^ and decays away from A. 
One can also write an Airy type asymptotic near |x| = A, but we do not need 
it here. The bounds of Lemma 15.11 follows from standard WKB bounds as in 
jl] and well-known formulas for the Hermite functions, see pQ. 

An immediate consequence of the above lemma is that all bounds in The- 
orem 0] are sharp in one dimension. In what follows we construct examples 
which show that the bounds in Theorem 0] are also sharp in higher dimension. 
For a positive integer N we consider eigenfunctions which correspond to the 
eigenvalue n + 2N. We define A > by 

A 2 = n + 2N 

5.1. Concentration in a tube in D . We set x' = (X2, ■ ■ ■ , x n ) and consider 
the eigenfunction 

v{x) = h N (xi)e~^ x '\ 
which is concentrated in the tube 

T = {|xi|<A, \x'\ < 1} 

We have 

\\v ||l 2 ~ 1 

and, if 1 < p < 00 

i_ 1 

IMlLP(mD^) ~ Af 2 . 
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5.2. Point concentration in D Q . We assume that N is even, otherwise the 
arguments need a small modification. We consider the set / of indices 

I = {a even; \a\ = N, «j > N/An) 

Then we consider the eigenfunction 

n 

a£l i=l 

which concentrates in the ball B(0, A^ 1 ). Since 
it follows that 

||u|| £a ~ A" -1 

On the other hand, by Lemma 15. 1\ for a 6 / we get 

h ai {xi) ~ A~5, \xi\ < A -1 

Summing up, we obtain 

h(x) ~ |J|A-f ~ A 2(n-l) A 2(n-l) ; % £ ^ ^1) 

therefore 

||/i||ip(fl(o,A-i)) ~ A-iA-^A 2 ^" 1 ) 

Thus 

A i-n( 3 --)|| v ||^^^ > i^ii^ 
which one should compare to (fTUj) . 

5.3. Point concentration in ZA,. We fix r with rei G DJ 11 and construct an 
eigenfunction which concentrates in B(rei, c2 J A _1 ). Let iVo be so that 

,»r A 2 -r 2 . 1 

A^ < 1. 

This implies that 

|JV | ~ 2- 2j A 2 

We consider the set of indices 

I = { a \ a' even; \a\ = N, a>i > N - N , a- > N/An} 
and its subset 

J = {ael; h ai (r)>±(\ 2 -r 2 )- 1 *} 

Lemma f5 . 1 1 shows that near r the functions h ai oscillate at a frequency of the 
order of 2~ J A. Then, changing r by no more that 2 J *A -1 , we can insure that 
\J\ ~ |/|. We define 



aGJ i=l 

Since 



|/|~iV 
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we have 

\\ V \\ L2 ~ |j|s ~ 2 -(™- 1 )iA n ~ 1 
On the other hand, using Lemma 15.11 and the definition of J we compute 



n-l 



v{x) ~ 25A"5iV * \J\ ~ A-f2Vi2- 2 ( n - 1 ) J 'A 2 ( n - 1 ), z £ 5(re 1; c2 J A~ 1 ) 
This gives 

therefore 

Al-™(l4) 2 ^-"(M)|K;|| LP(Drt) >K;|| i2 . 

5.4. Concentration on a tube in DJ 1 *. Here we construct eigenfunctions 
which concentrate in the tube 



T = {x£Df l - \x'\ <2-2} 



•tint 
3 

Now we consider the set of indices 



/ = { a - \a\ = N, a' even, |o/ - 2 j \ < c2 j } 
with small C. This has size 

|/| ~ 2 j(n_1) 

We want to consider a subset J of I of comparable size and set 

n 

The difficulty is that we want to avoid cancellations in this summation. For 
a £ I the functions h a i{x'^) are of size 2 _J//4 in an interval of length 2~ J / 2 
centered at zero. It remains to insure that the functions h ai {x\) have the 
same sign within most of the tube. For this we need to check that their phases 
are coherent, i.e. they differ essentially by a constant. 

Set /i = y/1 + ai. We must compare the functions s~(x) in the range 

A - \x\ ~ 2- 2j X, X 2 -fi 2 ~ 2 3 



For this we compute 



d d , . 2A 

- S -(x) = . ~ 2 J 



dji dx M vA 



Integrating this with respect to x and then with respect to fi we find that for 
x and /ii, /12 as above we have 

s m (X> - s~ 2 (x) = C(ni,fi 2 ) +0(2 _J A|/xi - /i 2 |) 

which yields 

l^iO) ~ s ^( x ) -C(fn,fi 2 )\ <c 
If c is small enough then the phases s~ (x) are close modulo constants, therefore 
we can choose a subset of indices fi of comparable cardinality so that the phases 

19 



are close modulo 2irZ. This leads to the subset J of I corresponding to this 
restricted set of a±s. 

Now we can complete the computation. On one hand 



\\v\\l2 ~ 2 2 3 
On the other hand, for most x G T we have 

71 + 3 ■ 

Since |T| ~ 2 z~ 3 \ this yields 

|M| LP(T) ~ 2-^A^-^'A-22-^ 

Hence 
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